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ABSTRACT 

Double peaked broad emission lines in active galactic nuclei are generally considered to 
be formed in an accretion disc. In this paper, we compute the profiles of reprocessing 
emission lines from a rclativistic, warped accretion disc around a black hole in order 
to explore the possibility that certain asymmetries in the double-peaked emission line 
profile which can not be explained by a circular Keplerian disc may be induced by 
disc warping. The disc warping also provides a solution for the energy budget in 
the emission line region because it increases the solid angle of the outer disc portion 
subtended to the inner portion of the disc. We adopted a parametrized disc geometry 
and a central point-like source of ionizing radiation to capture the main characteristics 
of the emission line profile from such discs. We find that the ratio between the blue 
and red peaks of the line profiles becoming less than unity can be naturally predicted 
by a twisted warped disc, and a third peak can be produced in some cases. We show 
that disc warping can reproduce the main features of multi-peaked line profiles of four 
active galactic nuclei from the Sloan Digital Sky Survey. 
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1 INTRODUCTION 

A small fraction of active galactic nuclei (AGN) 
show double-peaked broad emission line profiles 



pre 

jEracleous fc Halpernl 1 1994 l2003l : IStrateva etUT l2003l : 
iGezari. Halpern fc EracleousI 120071 ). The possibility has 
been considered for a long time that at least some of these 
hnes arise directly from the accretion discs assumed to 
feed the central supermassive black holes. The Hq profile 
observed in the spectrum of Ar p 102B was mode l led as 
arisi ng in a relativistic disc by IChen fc Halpernl (|l989l ) 
and IChen. Halpern fc Filippenkd (|l989l ) in an attempt to 
explain the asymmetric double-peaked profiles. In these 
stationary circular relativistic disc models, Doppler boost- 
ing will make the blue peak of the profile higher than the 
red one. However, iMiller fc PetersonI (|l990f ) questioned the 
relativistic disc explanation for Arp 102B, by showing that 
at least at some epochs the profile asymmetry is reversed, 
with the red peak higher than the blue one in contrary 
to the prediction of homogeneous, circular relativistic disc 
models. This type of profile asymmetry has been observed 
in some other double-peaked sources (jStrateva et al.l 
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l2003l : iGezari et al] l2007f ). Thus, other scenarios for the 
profile asymmetry were considered later on, including 
various physically plausible processes in the accretion flow, 
such as spiral shock s , ecc e ntric d iscs, bipolar outflows 
(IChakrabarti fc Wiital Il993l. Il994l: lEracleous et al.l Il995l; 



Storchi-Bergmann et al.l 119971 : IZheng. Sulentic fc Binettd 
199(3 : IVeilleux fc Zheng 199H). as well as a bina r y black 



hole system dB cgclma n. Blandford fc Reed Il980l : iGaskelll 
1 19961 : IZhang. Du ltzin-H acvan fc WaneJl2007bl V 



In addition to the line proflle asymmetry, the energetic 
budget for the line emission region has not been solved so 
far for the disc model. It is demonstrated that the local vis- 
cous heating in the accreti on disc is not sufficient to explain 
the o bserved Ha intensity (|Chen et al.ll 19891 : lEracleous et al.l 
[l997|), thus reprocessing of ultraviolet (UV) /X-ray contin- 
uum from the inner accretion disc is required. However, 
only a very small fraction of radiation from inner disc is 
expected to intercept the outer part of the disc in the case 
of a fiat geometrically thin disc. Two type of processes have 
been proposed to increase the fraction of l i ght in cident on 
to the disc emission-line region. IChen et al.l (| 19891 ') proposed 
that the inner part of the accretion disc is hot and be- 
comes geometrically thick due to insufficient radiative cool- 
ing, and the X-ray emission from such an ion-supported 
torus is responsible for such energy input. However, dou- 
ble peak emitters do not always have a low accretion rate 
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l|Zhang. Dultzin-Hacvan fc Wansj [200781 : iBian et al.ll2007h . 
Furthermore, a detailed calculation of advection dominated 
accretion flow (ADAF) models suggests that only a small 
porti on of X-rays actuall y hit the emission-line portion of the 
disc (ICao fc Wandbood ). On the other hand. ICao fc Wand 
(|2006l ') proposed that a slow moving jet can scatter a sub- 
stantial fraction of UV/X-ray light from the inner accretion 
disc back to the outer part of the disc. While this may be 
the likely case for radio-loud double-peaked emitters, the 
majority of double-peaked emitters are radio quiet. 

On the other hand, warped accretion discs are be- 
lieved to exist in a variety of astrophysical systems, in- 
cluding X-ray binaries, T Tauri stars, and AGN. The 
dynamics of warping accreti on discs has been discussed 
by a number of authors ([Papaloizou fc Pringld 19831: 
Papaloizou fc Lin!'l995!: 'O eavid1l999l : iNelson fc Papaloi^ 
1999, 2000; NaYakshin 200a). Four main mechanisms for 
exciting/maintaining warping in accretion discs have been 
proposed, including: tidall y induced warping by a compan- 
ion in a binary sy stem ( Terguem fc BertoutI 119931 . 119961 : 
iLarw ood ct al. 1996), radia tion driven or self-ind ucing warp- 
ing (LPringlc 1996, _1_997: Malonev. Begelman fc Pringlc 1996: 
Malonev fc BegelmanI 119971 : iMalonev. Begelman fc NowakI 



Malonev fe ijegelmanl lltJa/l : IMalonev. ijegelman fe INowaKl 
19981). magnetically driven disc warping (|Lail 1 19991 . l2003l : 



Pfeiffer fc Lail |2004|). and fram e dragging driven warp- 
teardccn fc PettersonI 1 19751: iKuniar fc Pringld 1 19851 : 



m^ : rr^'j- 

lArmitagc fc Nataraian 1 119991 ). As a result, the accretion 
disc in some AGN should be nonplanar. Evidence for the 
existence o f warped discs in AGN has been found by ob- 
servations (jHerrnstein et al.l l2005l : ICaproni et al.l 120071 ) . A 
warped disc makes it possible for the radiation from the 
inner disc to reach the outer parts, enhancing the repro- 
cessing emission lines. Thus, this warping can solve the 
long-standing energy-budget problem because the subtend- 
ing angle of the outer disc portion to the inner part has 
been increased by the warp. It can also provide the asym- 
metry required for the variations of the emission lines. The 
effect of disc warping on i ron line profiles treated rela tivis- 
tically have b een m a de bvlHartnoll fc BlackmanI (|200d ) and 
ICadez et all (|2003l ). iBachevI (|l999l ) studied the broad-line 
H/3 profiles from a warped disc. Nevertheless, his treatment 
is non-relativistic. A relativistic treatment of warped disc on 
the Balmer lines which can be compared with observations 
is still not available in the literature. 

In this paper, the broad Balmer emission due to re- 
processing of the central high-energy radiation by a warped 
accretion disc is investigated. Our main purpose here is to 
describe the effect of disc warping on the line profiles. A 
parametrized disc model is adopted. In !J5] we summarize 
the assumptions behind our model, and present the basic 
equations relevant to our problem. We present our results 
and the comparison with observations in ^ and conclusion 
and discussion in 34] 



2 ASSUMPTIONS AND METHOD OF 
CALCULATION 

In this paper, we will focus on how a warped disc affects 
the optical emission-line profiles, such as Ha and H/3. A 
geometrically thin disc and a central point-like source of 
ionizing radiation around a black hole are assumed. 



to observer 




Figure 1. The geometry used in this paper. The warped accretion 
disc is modeled as a collection of inclined rings, each defined by 
two Eulerian angles f3, 7 and the radial distance R. The normal to 
the equatorial plane is aligned with the Z direction. The direction 
to the observer is denoted by i that is the angle between the 
direction of line of sight and the Z-axis. 



2.1 Geometrical considerations 

The disc can be treated as being composed of a series of 
concentric circular rings of increasing radii, and laying in 
different planes. The rings interact with each other via vis- 
cous stresses. Each ring is defined by two Eulerian angles 
f3{R,t) and 'y{R,t) (see Fig.l) at radius R, where f3{R,t) is 
the tilt angle of the disc with respect to the normal to the 
equatorial plane, while the twist angle 7(-R, t) describes the 
orientation of the line of nodes with respect to a fixed axis 
in the equatorial plane. Note that _R is a purely radial co- 
ordinate and not the cylindrical one. At each radius R from 
the center, the disc has a unit tilt vector l{R,t) whic h varies 
contin uously with radius R and time t. Following IPringlg 
(|l996r ). in an inertial frame OXYZ centered on the disc, the 
vector l{R,t) is given by 

i = (cos 7 sin /3, sin 7 sin /3, cos /3). (1) 

We define the normalized vector towards the observer as 



iobs — (sin 1,0, cos i). 



(2) 



where i is the angle between the line of sight and the normal 
to the equatorial plane lie in the XZ plane. 

We define the coordinates on the surface of the disc as 
{R, (f)) with respect to a fixed Cartesian coordinate system 
{x,y,z), where </> is the azimuthal angle measured on the 
disc surface in the direction of flow, with cf) — 77/2 at the 
ascending nodelj The position vector of a point on the ring 
of radius _R is R = Rbr, where 



ck 



cos 7 cos j3 cos — sin 7 sin (; 

sin 7 cos P cos (j) + cos 7 sin i; 

— sin l3 cos (j) 



(3) 



is the radial unit vector. Note that I-br — 0, since points on 
the ring lie in the plane that passes through the origin and is 



This definition of (j> differs by 7r/2 from that of lPrinelel l ll996l) . 
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orthogonal to I. Now {R,4>) are (in general) non-orthogonal 
coordinates on the surface of the disc. The element of surface 
area is 



dS 



\dR 



AR 



dR 



A(j) 



(4) 



by calculating these partial derivatives directly and after 
some algebra which simplifies to 

dS = [Z + {R0 cos (f) + R'y' sin l3 sin (l>) sr] RdR d(f>, (5) 

where the primes indicate differentiation with respect to R, 
and hence, since I ■ en = 0, that 

|d5| = [l + {Rf3' cos (ji + R-/ sin p sin (PfY^'^ RdRd(l}. (6) 

The radiation flux L received from the central parts by the 
unit area is proportional to 



L oc 



|(_R/3 cos (fj + Ry sin 13 sin (j})\ 
R^ [l + (i?/?'cos<;/) + i?7'sin/3sin(/))^] 



211/2' 



(7) 



if the element in question is not obscured by other interior 
parts of the disc, otherwise the flux is set to zero (see be- 
low). The Balmer lines in the illuminated surface is formed 
as a result of photoionization by UV/X-ray radiation from 
inner accretion disc. We assume that Balmer line intensity is 
proportional to the radiation being intercepted by the disc, 
thus the line emissivity e on the disc surface can take the 
form: e oc L. A detailed treatment of line emission requires 
solving the vertical structure as well as the radiative trans- 
fer in the disc explicitly, which is beyond the scope of this 
paper. We also assume that the line emission is isotropic in 
the comoving frame. 

An alternative useful formalism being used to describe 
the shape of the warped disc is introducing the scale-height 
of the disc which is the displacement in the direction nor- 
mal to the equatorial plane. The scale-height denoted by 
/i(r, <p) is a function of radius r and azimuthal angle tp in 
cylindrical coordinates with respect to Cartesian coordinate 
system (X,Y,Z). The two formalisms described above can 
be related by the formula 



h{r,ip) — — _Rsin/3cos(^ 

= — _Rsin/3cos7/;/-\/l — sin^ /3sin^ ■;/) 
= —rta.nj3cos^. 

Where tp = tp — j and satisfies 
tan<^ = tan i/; cos /3. 



(8) 



(9) 



Since it is difficult to calculate analytically the disc dis- 
tortion induced by combined action of external torques men- 
tioned above, here we do not consider the disc dynamics, and 
merely parametrize the disc geometry to capture the main 
characteristics of emission-line profiles from warped discs. 
In this paper, a steady disc is adopted, and the form of the 
two angles 7, /3 as a function of radial R are given by 



7 
/3 



70 + nie 



713 sm 



1 ( JlZ 

2 V Rout 



R — Rin 



Ri, 



(10) 
(11) 



Where Rin and Rout are the disc inner and outer radii, 
Tii {i = 1,2,3) are three free parameters used to describe 
the warped disc geometry, 70 is a constant which is used 



to describe the longitude of the observer with respect to 
the coordinate system of the disc. This angle is mathemati- 
cally equivalent to changing observational time for possible 
precession. Thus, the line profiles with different 70 can be 
compared with the data observed at different times. 

2.2 Photon motions in the background metric 

The propagation of radiation from the disc around 

a Kerr black hole and the particle kinematics in 

the disc were studied by many authors (|Carteil Il968l : 
iBardeen. Press fc Teukolskvlll972l 'l. We review properties of 
the Kerr metric and formulae for its particle orbits, and 
summarize here the basic equations relevant to this paper. 
Throughout the paper we use units in which G = c = 1, 
where G is the gravitational constant, c is the speed of light. 
In Boyer-Lindquist coordinates, the Kerr metric is given by 



(12) 



ds^ 


= -e^^dt^ + e"^{d<j> - ujdtf + ^dr^ 


e^^ 


= EA/A, e"'* = sin^ SA/E, lv = 2Mar/A 


E = 


r^ + a cos^ 0, A = r^ -f a^ - 2Mr, 


A = 


(r^ + a'f -a'^sm'e. 



Here M, a are the black hole mass and specific angular mo- 
mentum, respectively. 

The general orbits of photons in the Kerr geome try can 
be ex pressed by a set of three constants of motion (jCarterJ 
Il96a ). Those are the energy at infinity _E, the axial com- 
ponent of angular momentum _EA, and carter's constant 
Q {=q E ). The 4-momentum of a geodesic has components 

Pp = (pt, Pr, P9, P0) = i-E, ±eVR/A, ±EVe, E\), (13) 
with 



R = 

e = 



A^ 



4 , / 2 

r + [a 



222 
q + a cos ( 



- q^) r^ + 2M [g^ + (A - af] r - 
A^ cot^ 9 . 



2 2 
a q 



From this, the equations of motion governing the orbital 
trajectory ca n be obtained. The motion in the r-8 plane is 
governed by (JBardeen et al.lll972l ) 



dr 



VW) 



de 



VW) 



(14) 



where r^ and 9e are the starting values of r and 6. The 
4> coordinate along t he trajectory is calculated by (jWilkinsI 
ll972l : IViergut3ll993l ) 



d<t) 



\de 



+ 



(2ar - Aa^)dr 



sin^e.jQ{e) Jr, A^/R{r) 
Consider the orbit equation in the form 
dr ^^f de ^ p 



(15) 



(16) 



where rma is the radius of the marginally stable orbit for 
a maximal Kerr black hole, a = 0.998. Then 9 can be ex- 
press ed in terms of P as follows (jCadez. Fanton fc Calvanil 
Il998h 



cos9{r) 



M+ 



cn(a^^^ + fi^_P ± K{m)\m) 



(17) 
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where m = 'a ^ 'a , fJ-± are defined by 



2 1 



and 



X' + q' 



a^)' + 4aV 



1/2 



T (A' + g - a 



A'(m) =cn"'(0|m). (19) 

Where K(m) is the complete elliptic integral of the first 
kind. 

For photons emitted from radii (>100rj) propagating 
to infinity the differences between Kerr and Schwarzschild 
black holes are indistinguishable. Neglecting terms of or- 
der 1/r^ and higher, the Kerr metric reduces to the 
Schwarzschild metric. The formulae for the particle orbits 
in this case are simplified. Here we summarize the ba- 
sic equations relevant to this paper. The line element in 
Schwarzschild geometry is written as follows: 

1 



ds^ = -(l-2/r)dr + - 



-Ar + rdd^ + rsm-dA4>' . (20) 



1-2/r 

The observer is assumed to be located at {ro,'&o,4'o)- 
For a approximately spherical background, each inclined 
ring of the disc can be considered lying in the equatorial 
plane. Owing to a spherically symmetric metric there is no 
favored direction for the black hole; in this sense, i9o can be 
acted as the angle between the observer and the normal to 
the disc and is determined by taking the dot-products of the 
unit vector I of the disc with the normalized vector iobs to 
the observer. By definition, we get 

cos i?o = cos 7 sin /3 sin i -f cos /3 cos i (21) 

The angle 00 can be obtained by setting the mixed 
triple product of three vectors iobs,i and ea equal zero, 
(fobs I ea) = 0. From this equation, the angle is determined 
by 



tan( 



sm 7 sm i 



I2OOII ) 



• (22) 

sin P cos i — cos /3 cos 7 sin i 

The photon motion equations are given by (|Lu fc Yul 



dr 



y/Rir 



d4>- 



= ± 



di9 



Adi9 



(23) 



(24) 



where re, i9e and (j>e are the starting values of r, i9 and (f). 

Define fi = cosi9, then in equations (|23|) and (|24|) the 
integral over ^ can be worked out in terms of a trigonometric 
integral 



d& 

./2\/eM 



\M 



dfi 



^g2 - (A2 + q^) fi^ 



V^+q'- 



:sin ^ (m/m-i-), 



(25) 



^/2 sin^^yeo?) 



Xdfi 



= ±sin 



{i~^^^)./eM 



^^(i-M^)- 



(26) 



Where ^ ^ < ^+ , /i+ = \nf7i^^ + <f) and 0^(At) = 
q — (A +q )n . The integral over r can be worked out with 
inverse Jacobian e lliptic integrals (see e.g. ICadez et al.lll998l : 
( Win fc Wa^l2007l '). 

Although our calculation is limited to Schwarzschild 
metric for photons travelling to infinity, for the radii (>100 
Tg) taken into account in this paper, the results are essen- 
tially the same for the Kerr black hole. On the other hand, 
for the light from inner part of the disc, a maximal Kerr 
metric is still used. 



2.3 The observed line flux 

Due to relativistic effects, the photon frequency will shift 
from the emitted frequency Us to the observed one Uo re- 
ceived by a rest observer with the hole at infinity. We in- 
troduce a g factor to describe the shift which is the ratio of 
observed frequency to emitted one: 

g = !/o/!^c = (p ■ Uo)/{p ■ Uc) 

= (l-nA)-^^!^, (27) 

where Q, = r^ ' is the Keplerian angular velocity, p, Uo, Ue 
are the 4-momentum of the photon, the 4-velocity of the 
observer and the emitter, respectively. 

The specific flux density -Fo(i^o) at frequency Vo as ob- 
served by an observer at infinity is defined as the sum of 
the observed specific intensities Io{i^o) f rom all p arts of the 
accretion disc surface, which is given by (|Cunninghamiil975l ) 



Foil^o) 



Io{lyo)dflobs 



g Ic{uc)dQob 



(28) 



where dilahs is the element of the solid angle subtended by 
the image of the disc on the observer's sky and we have made 
use of the relativistic invariance of lu/iy^, v is the photon 
frequency measured by any local observer on the path. 

Icii'a) is the specific intensity measured by an observer 
corotating with the disc, and can be approximated by a 5- 
function, Ie{v'e) = e5{i'e — Ve), where e is the emissivity per 
unit surface area. From well-known transformation proper- 
ties of 5-functions we have 5{u'c — Vc) = g5{vo — gve), using 
this in equation (|28|) . we obtain 



(29) 



In order to calculate the integration over df^obs, using 
two impact parameters a and /?, measured relative to the 
dir ection to the centre of the blac k hole, firstly introduced 
by ICunningham fc BardeenI (| 19731 1 is convenient. They are 
related to two constants of motion A and q by equations 

Q = -A/sim?o, /3 = ± (g^ - A^ cot^ i9o) ^''^ (30) 

The element of solid angle seen by the observer is then 



Foifo) = / eg*5{iyo - gUa)dVLoh 



dfiobs ~ 



dadj3 



d{X,q) 



dr dg, 



(31) 



ro/3sini9o d{r,g) 

where ro is the distance from the observer to the black hole. 
Substituting equation (|31|l into equation (|29[l gives the 
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desired result: 



Fo{v^) = 



nate form of equation (|27 



3/3sin j?o 



eg'^Sivo - gvc) „, ' . dr dg. (32) 



In the calculation of the total flux of reprocessing emis- 
sion line, shadowing of the elements by the inner parts of the 
disc must be taken into account. This has been done as fol- 
lows. The contribution of each element denoted by (r, ip) to 
the total line flux from the disc is calculated by equation (|32p 
only in case that there does not exist another element with 
cos 6i > cos 9{ri, tpi) tracing the trajectory of photon, and is 
zero in the opposite case. Here cos^i = h{ri < r,ifii)/ri, tpi 
and cos6{ri,ipi) are calculated by equations (|15|l and (|17p . 
respectively: 



ifi ^ ip- 



2ardr 



, ipi) = ^+cn(a-\/^^ + M^P — K{m)\'m) 



(33) 



(34) 



Equation (I33|l is obtained by setting A = in equation (|15p 
as a approximation for a pointlike source assumption. In the 
calculation, the disc is divided into 1000 rings logrithmically 
spaced in radial direction , and we checked every increment 
of r for the inequality along the ray path. In the calculation 
we neglect the contribution from the lower surface of the 
disc or the higher order images of the emitter due to the 
fact that at low observer inclination angles the extra flux is 
small compared to the direct images. And a more correct 
calculation will need to handle these effects correctly and 
ma y give up to a f actor of two change in central region (see 
e.g. lViergut3ll993l ). 



2.4 Method of calculation 

With all of the preparation described in the previous section, 
we now turn to how to calculate the line profiles numerically. 
We divide the disc into a number of arbitrarily narrow rings, 
and emission from each ring is calculated by considering its 
axisymmetry. We shall denote by ri the radius of each such 
emitting ring. For each ring there is a family of null geodesies 
along which radiation flows to a distant observer at polar 
angle i?o from the disc's axis. As far as a specific emission 
line is concerned, for a given observed frequency fo the null 
geodesies in this family can be picked out if they exist. So, 
the weighted contribution of this ring to the line flux can 
be determined. The total observed flux can be obtained by 
summing over all emitting rings. 

The main numerical procedures for computing the line 
profiles are as follows: 

(i) Specify the relevant disc system parameters: 
Tin, rout, ni, 712, ns and i, 70 . 

(ii) The flux from the disc surface is integrated using 
Gauss-Lege ndre in tegration with an algorithm due to Ry- 
bicki G. B. (|Press c t al. 199^). The routine provides abscis- 
sas Ti and weights lji for the integration. 

(iii) For a given couple (ri, g) of a ring, the two constants 
of motion A and q are determined if they exist. This is done 
in the following way. The value of A is obtained by an alter- 




(35) 



the value of q is determined for solving photon trajectory 
equation (|23[) . Then the contribution of this ring on the flux 
for given frequency Uo with respect to g is estimated. Varying 
g, this step is repeated, thus the fiux contribution of this ring 
to all the possible frequency is obtained. 

(iv) For each g, the integration over r of equation (|32[) can 
be replaced by a sum over all the emitting rings weighted 
by a;i 



Fo(;/o) 



4 
qeUo 



d{X,q) 



i— 1 



dir,g) 



(36) 



The Jacobian [d{X, q)/d{r, g)] in the above formula was eval- 
uated by a finite difference scheme. In the calculation the 
self-obscuration along the line of sight are also taken into 
account. This has been done as follows. The contribution of 
each element denoted by (n, 5) to the total line fiux from the 
disc is calculated by equation (|36|) only in the case that there 
does not exist another element with cos6'i > cos6{ri,ipi) 
along the ray path, and is zero in the opposite case. Here 
cos Si = h{ri,ipi)/ri, cos9{ri,ifii) calculated by: 



cos e{ri,(pi) = 



fi+ sin(7r - KM - ^/A^-Fq^P), 



M+ 



sin(KM±yA2+q2P), 



(37) 



where the first equation describes the case with one turning 
point in 9 component along the geodesic, tpi can be calcu- 
lated by equation (|24[l and equation (|26p. and P,KM, fi+ 
satisfy 



dr 



VW) 



de 



VW) 



rr, T . _1/COSZ, 

KM = sm ^( ) 



P+ 



VWi^^+^ 



From the above formula, one determines the line flux at 
an arbitrary frequency Vo from the disc. The observed line 
profile as a function of frequency Uo is finally obtained in 
this way. 



3 RESULTS 

In the accretion disc model, the double-peaked emission lines 
are radiated from the disc region between around several 
hundred gravitational radii to more than 2000rg; here r^ 
is the gravitational radius and the widths of double-peaked 
lines range from s everal thousand to nearly 40,000 kms~^ 
I Wang et al.ll2005l ). In our model, all of the parameters of 
the warped disc are set to be free. The reasonable range are: 
disc radius from lOOrj, to 2000r-g, ni from to Stt, 712 from 2 
to 4, 712 from to 30° . This model has the same number of 
free-p arameters as an eccentric disc model (jEracleous et al.l 
[1995I). In the range of frequency from 4.32 to 4.78 in units 
of 10^'^Hz, 180 bins are used. Considering the brodening 
due to electron scattering or turbulence, all our results are 



6 Wu, Wang and Dong 



Wt 



-i4H 



C" 



1 1 1 m i nHi i M >■ ■ ■'! 



^ 



€ 



n I I I iF- mi ii i>F n ii 



3 



h n i'^ii I il iF r >■ » < II I il f FF rm ii 



II II i K ii'lii I i| I' r Fi j i II I i| f FF F S M iil II I 'r ' 



^ 



nu ll i i T . M i "ti i i r t 



«> 



O 



€^ 



~iDO>-aD a 3i» im) 



-JMO^SK fl sea locc 



-lOBJ^HM a «0 10(H 



«. 



W 









ih jr i i-^-, - »_ . ^ h 




11:^ 




k f I I nil I I 




J^-i— — i—l^A- 



<a» ISM ann eaa esm <s» u» sm S3» esK ATm- u» u» laa exa ttst u» asn tna axa em g&» imu tsa 



Wavelength (A) 

Figure 2. The images (upper panel) of the illuminated area of the disc and Ha line profiles (lower panel) computed by our code for 
twist warped disc case for i = 30°, ni = 2.l7r, n2 = 3.0, n^ = 10°, the disc zone is from Hin = lOOrg to -Rout = lOOOrg. The profiles 
from both prograde disc (black line) and retrograde disc (dotted line) are shown. The longitude of the observer 70 with respect to the 
coordinate system of the disc takes steps of 30° from 0° to 330° (from top left-hand side to bottom right-hand side). All the results are 
smoothed by convolution with 3a Gaussian (corresponding to 500 kms^^). 



ii» 



sxt 



-4«4 



r 



I t I nH i L Fii iFi H f h i ; 



>3 



C 



i. lFF lFii m II ifFIFli 



5> 



V 1(9 



11 I I I > F I ill 1 I I I I I F . Fl i t I I I I I F F F Fl m I I I I ■ ■ t 

<y\cy 



V9 



m ill FU FhF m lll ••> 



o 



^ 



-iDso-oai a SCO IMD 



-!flM^3l» B BOO lOCC 



^HHn-^aH II ;cD -taoa 



J9 



* 1 1 . i ti 





. i . . Ti I , . . i > I I 




■ ■ >■ ■ ' .^ > ■ - ■ > > ■ - ■ >. 





..ill. N I ' li 




tsK' tstt CTDC eaa ssa tsw ime e«» ema esu ATm mod exu laa exu oca uu asas it^ t»xi £iiso nna imK ts>d4 



Wavelength (A) 



Fieure 3. As in FiEl2lbut with i = 15°. n.i = 20° 



Emission from warped discs 7 




6300 6500 6700 6900 6300 6500 6700 6900 6300 6500 6700 6900 6300 6500 6700 6900 6300 6500 6700 6900 6300 6500 6700 6900 

wavelength 

Figure 4. Comparison of the Hq line profiles generated by our code for warped disc cases with different phase amplitude (described by 
parameter ni), ni = 0, 7r, 27r, Stt (from top to bottom panels). The other parameters are: i = 30°, 712 = 3.0, na = 10°, Rin = lOOrg 
and -Rout = 1200rg. The longitude of the observer 70 with respect to the coordinate system of the disc varies from 0° to 300° with a 
step of 60° (from left- to right-hand side). 



smoothed by convolution with Scr Gaussian (corresponding 
to 500 kms^M. 



3.1 Line profiles from twisted warping discs 

We ex tend the numerical code developed by IWu fc Wang 
((20071) to deal with the warped discs. The images of the illu- 
minated area and the Ha line profiles computed by our code 
for a steady-state twisted warping disc are shown in Figs [5] 
andO The disc zone is from Rin = lOOrg to 7?out = lOOOrg, 
the other parameters are i = 30°, ni = 2.l7r, 712 = 
3.0, ns = 10°, and i = 15°, rii = 2.l7r, 712 = 3.0, n-i = 20°, 
respectively. The longitude of the observer 70 with respect 
to the coordinate system of the disc from top left-hand side 
to bottom right-hand side varies from 0° to 330° with a 
step of 30° . The image of the illuminated portion of the disc 
has a shape similar to a one-armed disc. The image area of 
the disc being projected on the observer's sky from the far 
side of the disc due to a little angle between the tilt vector 
1{R, t) and the line of sight significantly larger than the near 
parts are indicated in Figs [5] and O In our code, one can 
specify an arbitrary grid for r or 0. In our calculations, all 
the images contain 400 x 360 pixels in this paper. The line 
profiles drawn in black are restricted to prograde disc mod- 
els (relative to the normal to equatorial plane or black hole 
spin), the effects of a retrograde disc (identical to the pro- 
grade case with a reversed spin) on the line profiles should 
also be taken into account if the orientations of the AGN 
are isotropically distributed in space, those are drawn with 



dotted line in the plots. The profiles from a prograde disc 
and its corresponding retrograde disc is close to but not ex- 
actly reflective symmetric about v — because of relativistic 
boosting, the gravitational redshift and the asymmetry of 
the disc. Note also that there are triplc-pcakcd profiles simi- 
lar to those found by observations ( Vcillcux & Zheng 199l|) 
and the obvious variation in the red/blue peak positions. 



The influence of the phase amplitude (described by pa- 
rameter rii) on the line profiles for ni = (twist-free), 
TV, 2-K, 37r are shown in Fig. [4] The disc zone is from 



Rin 



1200rg, the other parameters are 



i = 30°, 712 = 3.0, 713 — 10°. The longitude of the ob- 
server 70 with respect to the coordinate system of the disc 
varies from 0° to 300° with a step of 60° (from left- to right- 
hand side). For twist-free and low phase amplitude discs, the 
asymmetric as well as frequency-shifted single-peaked line 
profiles expected in most cases are produced. For tii > 27r, 
there is a fairly large possibility for a triple-peaked line pro- 
file. The influence of the phase amplitude on the line profiles 
is remarkable. 



The fraction of light incident to the disc from the in- 
ner region is a function of 77i , 713 . We calculate the fraction 
received by a warped disc, the results are shown in Fig. [S] 
The fraction reaches 27.3 per cent for a twisted warp disc 
m = 2.07r, 712 =3.0, 713 =20°. 
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Figure 6. comparison of the Ha line profiles computed 
by our code (red line) with the observation of SDSS 
J084205.57+075925.5. The parameters are: i = 13°, nj = 
S.Ott, n2 = 3.0, n3 = 28°, 70 = 95°, the disc zone is from 
iJin = 250rg to iJout = lOOOr'g for a prograde disc. 



3.2 Comparison with the observations 

As the first step, our main purpose in this paper is to demon- 
strate the efTect of disc warping on the line profiles, while 
the detailed fitting of the observational data will be pre- 
sented in a future paper. From the sample of AGN with 
double-peaked Balmer lines drawn from the Sloan Digital 
Sky Survey (SDSS) spectroscopic data set (Shan et al., in 
preparation), we pick up four sources with peculiar profiles: 
SDSS J084205.57+075925.5, SDSS J232721.96-f^l52437.3, 
SDSS J094321.97-h042412.0, SDSS J093653.84-f533126.8. 
Two of them have explicit triple-peaked profiles, the oth- 
ers show a net shift of the emission lines to the red. These 
features are hard to interprete via a homogeneous, circu- 
lar relativistic disc model. The SDSS spectrum of one ob- 
ject (SDSS J232721.96-I-152437.3) is dominated by starlight 
from the host galaxy; we subtract the starlight and the 
nuclear continuum with the method as described in detail 
in Zhou et al. (2006); The other three spectra have little 
starlight contamination, the nuclear continuum and the Fe II 
emission multiplets are modelled as described in detail in 
Dong et al. (2008). The Fell emission lines are apparent in 
J0936-F5331. 

The observed broad line spectrum of sources SDSS 
J084205.57-f075925.5, 

SDSS J232721. 96+152437.3 both have a third peak. The 
line profiles (red line) calculated by our code compared with 
observations are shown in Figs [6] and [7] The disc zone is 
from i?in — 250rg to i?out = lOOOrg, and 712 — 3.0 for a pro- 
grade disc case. The other parameters are: i — 13°, ni = 
3.07r, ng = 28°, 70 = 95° for SDSS J084205.57-f 075925.5, 
i = 10°, ni = 2.67r, m = 22°, 70 = 88° for SDSS 
J232721. 96-1-152437.3, respectively. We have taken the ap- 
proach that the disk-like component should be fitted by eye 
to as much of the peaks match as possible. While our model 
does not reproduce the extended wing, the peak positions 
and heights are reasonably in agreement with the observed 
profiles. 

A retrograde disc is needed to match the line pro- 
files with the observations for SDSS J094321.97-f 042412.0, 



SDSS J232721. 96 + 152437. 3 
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Figure 7. comparison of the Ha line profiles computed 
by our code (red line) with the observation of SDSS 
J232721. 96-1-152437.3. The parameters are: i = 10°, ni = 
2.67r, 712 = 3.0, n-j = 22°, 70 = 88°, the disc zone is from 
-Rin = 250rg to -Rout = lOOOrg for a prograde disc. 



J093653.84-f533126.8. Fig. |8] shows the comparison be- 
tween the line profiles and the observed profile for SDSS 
J094321. 97+042412.0 with i = 17°, m = 2.l7r, 713 = 
24°, 70 = 202°. And the parameters in Fig. H for SDSS 
J093653.84+533126.8 are: i = 8°, m = 2.l7r, ns = 6°, 70 = 
117°. The disc zone is still taken from Rin = 250rg to 
^out = lOOOrg and 712 — 3.0. We see that the overall 
match except that in the two wings between the line profiles 
and the observations is good qualitatively. The deviation in 
two wings may, mostly due to the central point-like source 
adopted, diminish the contribution from the relatively inner 
part of the reprocessing region. 

The line profiles with respect to different longitude 
of the observer are corresponding to the observations in 
different times in reality. Thus the variation of the line 
profiles with respect to different longitude can be com- 
pared with successive observations. The images of the 
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Figure 8. comparison of the Ha line profiles computed 
by our code (red line) with the observation of SDSS 
J094321. 97+042412.0. The parameters are: i = 17°, ni = 
2.1-77, 712 = 3.0, na = 24°, 70 = 202°, the disc zone is from 
iiin = 250rg to iJout = lOOOr'g for a retrograde disc. 
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Figure 9. comparison of the Ha line profiles computed 
by our code (red line) with the observation of SDSS 
J093653.84+533126.8. The parameters are: i = 8° , rii = 
2.l7r, 712 = 3.0, na = 6°, 70 = 117°, the disc zone is from 
Rin = 250rg to -Rout = lOOOrg for a retrograde disc. 



illuminated area and the Ha line profiles for sources 
SDSS J232721.96+152437.3, SDSS J093653.84+533126.8 
are shown in Figs [lO] and 1111 In Fig. [lO] the parameters of 
the disc are: i = 10°, rii = 2.67r, 712 — 3.0, 713 — 22°, 
-Rin = 250rg, Rout = lOOOrg for a prograde disc. The longi- 
tude of the observer 70 from top left-hand side to bottom 
right-hand side takes steps of 10° from 30° to 140°. A ret- 
rograde disc with i — 8°, rii = 2.l7r, 722 = 3.0, 713 — 6°, 
-Rin = 250rg and -Rout ~ lOOOrg, the longitude of the ob- 
server 70 from top left-hand side to bottom right-hand side 
takes steps of 10° from 50° to 160° (shown in Fig. Ill|l . Long 
term monitoring of these sources should provide critical tests 
for the disc warping scenarios. 



4 SUMMARY 

We compute the reprocessing emission-line profiles from a 
warped relativistic accretion disc around a central black 
hole by including all relativistic effects. A parametrized disc 
model is used to obtain an insight into the impact of disc 
warping on the double-peaked Balmer emission-line profiles. 
For simplicity, we assumed that the disc is illuminated by a 
point-like central source, which is a good approximation, the 
line emissivity is proportional to the continuum light inter- 
cepted by the accretion disc, and line emission is isotropic. 
We find the following. 

(i) For twist-free or low phase amplitude (described by 
ni) warped disc, the asymmetrical and frequency-shifted 
single-peaked line profiles as expected are produced in most 
cases. The rarity of such sources suggests that a warping 
disc is usually twisted. 

(ii) The flux ratio of the blue peak to red one becoming 
less than unity can be predicted by a twisted warp disc with 
high phase amplitude. The phase amplitude has a significant 
influence on the line profiles. 

(iii) A third peak and the variation in blue/red 
peak position which have been found by observations 
iveilleux fc Zheng |[l99ll : IStrateva et al.ll2003l ') are produced 
in our model. 

(iv) The influence from retrograde disc on line proflles has 
a inverse effect compared with a prograde one. 

(v) The fraction of the radiation incident to the outer disc 
from the inner part can be enhanced by disc warping. The 
results are shown in Fig. \5\ 

We showed that warped disc model is flexible enough 
to reproduce a variety of line proflles including triplet- 
peaked line proflles, or double peaked proflles with addi- 
tional plateaus in the line wing as observed in the SDSS 
spectra of some AGN, while the models have the same num- 
ber of free parameters as eccentric disc models. While we 
leave the detailed fit to a future paper, the essential char- 
acteristics in the line proflles is in good agreement with ob- 
served line proflles. Future monitoring of line proflle vari- 
ability in these triplet sources can provide a critical test of 
the warped disc model. 
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